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Quantum manifestations of the dynamics around resonant tori in perturbed Hamiltonian sys- 
tems, dictated by the Poincare-Birkhoff theorem, are shown to exist. They are embedded in the 
interactions involving states which differ in a number of quanta equal to the order of the classical 
resonance. Moreover, the associated classical phase space structures are mimicked in the quasiprob- 
ability density functions and their zeros. 

PACS numbers: 05.45.Mt, 03.65.Sq 



The quantum manifestations of classical chaos, or 
quantum chaos, have received much attention in the re- 
cent past [l]. Besides random matrix theory, which de- 
scribes universal properties, Gutzwiller's trace formula 
describes the dual relationship between coherent sums 
over periodic orbits with sums over quantized states [H. 
The appearence of scars of periodic orbits in individual 
eigenfunctions and the associated scar theory are results 
that deserve special mention |3H5]. These results are by 
now routinely applied in a variety of important technolog- 
ical applications [6]. At the other extreme, for integrable 
systems, this dual relationship is more precise and collec- 
tive sums over the degenerate periodic orbits on rational 
tori give rise to individual states on quantized tori j?]. 

For systems close to integrable, the classical behav- 
ior is completely understood in terms of the celebrated 
Kolmogorov-Arnold-Moser (KAM) and the Poincare- 
Birkhoff (PB) theorems ji, The first one deals with 
quasiperiodic motion and the persistence of sufficiently 
irrational tori under perturbations. The second considers 
the fate of the (unperturbed) resonant tori, from which 
an even number of POs survive. In their vicinity, chains 
of islands of regularity surrounded by a chaotic separa- 
trix organized by a homoclinic tangle are formed. This 
structure and the associated transport are well described 
by the universal pendulum model of Chirikov [lo|, be- 
ing the (slow) Arnold diffusion [8] the controlling mech- 
anism of transport. Classical resonances are known to 
control many relevant processes, such as intramolecular 
vibrational relaxation or directional laser emission 
in nanooptics [l, [l2|. On the other hand, the analysis 



of the quantum counterpart has produced much fewer 
results [13-15]. 

In this Letter, we explore the quantum implications of 
the classical PB theorem. We find that they exist and 
can be unveiled by studying the mechanism by which 
two quantized tori of the unperturbed system interact 
to form the beginning of an island chain on a resonant 
rational torus. The subtle mechanism by which the quan- 
tum numbers of the unperturbed tori are exchanged in 
a Landau-Zener (LZ) pj| quasi-crossing is illustrate d by 
the exchange of zeroes of the stellar representation [17|. 
The splittings of a series of avoided crossings (AGs) ruled 
by a PB resonance in the correlation diagram is explained 
semiclassically. 

The model that we have chosen to study is the Harper 
map in the unit square, 

Pn+i = Pn^ksm{27rqn) (modi), 

<?n+i = - ^sin(27rpn+i) (modi), (1) 

where /c is a parameter measuring the strength of the 
perturbation. The quantized version is provided by the 
unitary time-evolution operator [Til, Il9| 

U = exp[iNkcos{27Tq)] exp[i7Vfe cos(27rp)], (2) 

with TV = {27rh)~^. The simplicity of this model allows 
extremely detailed calculations. The classical dynamics 
of this map follows the KAM and PB pattern as a func- 
tion of /c, depicted in Fig. [Ha)-(d). For small values 
k = 0.1 [panel (a)] all trajectories look regular at the 
scale of the figure, being confined to invariant tori. For 
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FIG. 1: (color online) Classical phase space for the Harper 
map dU) for k: (a) 0.1, (b) 0.155, (c) 0.2, and (d) 0.25. 
(e) Correlation diagram for the eigenphases of the quantized 
Harper map ([2|) with A/" = 60 as a function of the perturbation 
parameter k. Two series of avoided crossings of states with 
quantum numbers differing in 10 (red squares) and 6 (blue 
circles) units, respectively, are shown. 



k = 0.155 [panel (b)] two resonances of order 1:14 and 
1:10 around the center become important and the corre- 
sponding chains of islands are visible. Also, the separa- 
trix between central tori and those at corners shows as an 
appreciable band of stochasticity. As k reaches 0.2 [panel 

(c) ] the classical phase space continues evolving; the 1:14 
resonance disappear, swallowed by chaos in the main sep- 
aratrix, the 1:10 resonance grows in size, and a new 1:6 
resonance becomes apparent. Finally, for k = 0.25 [panel 

(d) ] only the 1:6 resonance survives. 

The quantum map can be diagonahzed resulting in a 
set of TV eigenphases, (j)n G [— 7r,7r), and eigenfunctions, 
\(j)n)- Plotted as a function of k the eigenphases produce 
the correlation diagram in Fig. [U^e). To interpret this 
diagram we rewrite the map using the Baker- Hausdorff 
expansion for non-commuting exponentials ^ 



21| as 



U = 



exp exp 



(3) 



where = Hq ^ kHi ^ ■ ■ ■ k^'-^Hn-i and V^''\k) is 
another operator with a power series expansion starting 
with the power k'^^^ . All terms of both expansions can be 



computed in terms of iterated commutators of cos(27r^) 
and cos(27rp). If commutators are replaced by Poisson 
brackets in the usual fashion, the equivalent classical ex- 
pansion is obtained [2l|. As an example the classical 



expansion to second order yields 



cos(27rg) + cos(27rp) 
2^r 



■ sin(27rg) sin(27rj9) 
(4) 



V^'^\p^q) — [^^^(^'^^) ^^^^(^'^^) + cos(27rj9) sin^(27rg) 

Further terms can be (laboriously) calculated. The im- 
portance of this rewriting for our purposes is that for 
moderate values of k the exact eigenphases can be in- 
terpreted as the levels of the integrable operator H^{k) 
with definite quantum numbers n. When multiplied by 
kN and wrapped around the unit circle, the wrapping 
produces AGs of states with different quantum number 
which become AGs of LZ type [l6| when the small terms 
of order in V are taken into account. This pro- 

cedure becomes eventually divergent reflecting the fact 
that the chaotic behavior of the map cannot be captured 
by any integrable approximation. However the proce- 
dure is well suited to describe the beginning of the PB 
process when a rational torus disintegrates in an island 
chain. This situation is clearly visible in the correla- 
tion diagram of the eigenphases, calculated for N = 60. 
There are 60 eigenlines all starting at zero for /c = 0, and 
then they spread in an fan-like fashion, first linearly with 
k. The slopes are given initially by the eigenvalues of 
J^o = [cos(27rg) + cos(27rj9)]/27r and they are ordered in 
increasing value of their quantum number. Half of these 
lines have positive slope corresponding to states localized 
around (p, q) = (0, 0) while the ones with negative slopes 
are centered on (p, g) = (1/2,1/2). We will restrict our 
analysis to the island formation in this latter case, the 
other one being equivalent since connected by a simple 
symmetry operation. 

Two interesting features arise from the analysis of 
the correlation diagram. First, the AGs always appear 
grouped in series or families, where the involved states 
differ by the same number of quanta. An. Two such 
families of AGs, corresponding to An = 10 and An = 6 
have been identified and marked with red squares and 
blue circles, respectively, in the correlation diagram of 
Fig. [H The detailed analysis of one of these AG is shown 
in Fig. [21 The AG involves states with niower = 4 and 
^upper = 10 in the An = 6 family and is a typical sharp 
two-level LZ AG [lej. The left panel of Fig.[2]is a blow-up 
ofFig.fflintherange A: = 0.185-0.21. Points (a)-(f), (b)- 
(g), (c)-(h)(d)-(i) and (e) and (j) show different stages of 
the AG, the minimum A^ being attained at point (d)-(i). 

More interesting is the evolution of the associated 
Husimi functions, which is shown, using the same la- 
belling, in the right panels of the figure. Recall that the 
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FIG. 2: (Left) Detail of the correlation diagram of Fig. [T] corresponding to the fifth element of the family of avoided crossings 
with An = 6 which involves states with quantum numbers ni — 4 and n2 — 10. 

(Right) Evolution of the corresponding quasiprobability densities given by the Husimi functions at selected points [(a)- (j)] in 
the diagram. Full circles mark the position of the zeros of the Husimi functions. 



Husimi distribution of a quantized torus is localized on it, 
and its quantum number is given by the number of zeroes 
enclosed by the torus [HI- The rest of the zeroes, in 
this case to a total of A/" = 60, is scattered in the comple- 
mentary phase space region. Far from the AC, at points 
[(a)- (f), (b)-(g)] the two tori maintain their individual- 
ity, but the zeroes start moving in a very characteristic 
way: the lower state n = 4 attracts 6 zeroes while 6 ze- 
roes start migrating outwards in the upper n = 10 state. 
The process culminates in panels (d)-(i) where the two 
states are now precisely localized on the rational torus 
but with maxima either on the stable or the unstable pe- 
riodic points given by the PB theorem (the upper state 
sits on top of the unstable points and the lower on the 
stable points). Similarly, the position of the zeroes alter- 
nate between the fixed points, localizing in the immediate 
neighborhoods. Beyond the AC, the original topology is 
restored, aside for the compulsory exchange of quantum 
numbers. We have confirmed this mechanism for all other 
ACs of this family and also for the 1:10 resonance. We 
can then conclude that this is a general mechanism by 
which two quantized tori interact via the periodic orbits 
- stable and unstable - that, according to the PB theo- 
rem, remain on the resonant torus. Even more interesting 
is the fact that the number of zeroes exchanged between 
the tori is precisely the order of the resonance. 

Let us now make the above arguments quantitative. At 
a LZ AC the splitting of the levels is given by the inter- 
action matrix element between the two involved states. 
Quantum mechanically it could be calculated from the 
eigenstates of the integrable hamiltonian i^(^) perturbed 
by V^^"^ in Eq. (|3]) and computed up to a given power of k. 
However, we are interested in a purely semiclassical cal- 
culation. To avoid the very cumbersome transformation 
to action- angle variables, in which the calculation would 
be trivial, we resort to a method developed in Ref. [22 1 



in connection with the calculation of tunneling rates in 
this same model. The resulting expression for the eigen- 
phase difference at the AC due to a 1 : r resonance is 
then obtained in terms of purely classical information as 



nA(t) = 



1 



exp(-zr(9) SIi:r{0)dO, 



Sli:r{0) being given by 



(5) 



(6) 



where /^~^^(/i:^, 0) symbolizes the action variable that is 
obtained by applying the inverse Poincare map to (/, 0) 
or, alternatively, the backward propagation with H from 
time t = r to t = 0. The way the calculations are done 
in practice is the following. First, the resonant periodic 
torus is located, and its action computed. This is done 
by using an analytical third order approximation to the 
Harper Hamiltonian in Eq. (jSj , followed by numerical re- 
finement propagating trajectories near the torus. Once 
the resonant torus is found, the corresponding angle vari- 
able is computed for a large sample of points, which are 
then back-propagated with the exact inverse map. (This 
is a trivial process since it just implies exchanging the p 
and q variables.) The associated action, corresponding 
to each of these points, is then calculated by numeri- 
cal propagation of the trajectories in a complete cycle. 
This provides the values of the perturbed action, from 
which A(j) is evaluated with the aid of Eq. (|5]). The cor- 
responding results for the 1:6 resonance are shown as a 
continuous line in Fig. [3l As can be seen, the agreement 
of this semiclassical calculation with the data (triangles) 
extracted from the correlation diagram of Fig. [T] is ex- 
tremely good, this giving full support to our previous 
qualitative arguments on the existence of a clear quan- 
tum manifestation of the classical PB theorem. 
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FIG. 3: (color online) Numerical eigenphase differences, A(f), 
for avoided crossings of the An = 6 family extracted from 
the correlation diagram of Fig. [T]^blue circles) as a function of 
the perturbation parameter for N: (gray A) 30, (blue 0) 60, 
(red x) 100, and (full line) semiclassical calculations based on 
Eq. Q for the 1:6 resonance. 
Inset: Same for the 1:10 resonance and N = 100. 



rationalize classical Hamiltonian chaos, exist. They can 
be unveiled by analyzing the interactions, or AGs split- 
ting in a correlation diagram, involving pairs of states 
which differ in a number of quanta equal to the order 
of a given classical resonance. This result has remark- 
able effects, since as a (perturbational) parameter varies, 
and a quantized torus passes through the resonance it 
looses (or gains) zeroes in its Husimi representation. At 
resonance the zeroes are located at the position of the 
periodic points predicted by the PB theorem, and the 
number of zeroes, lost or gained, is equal to the order of 
the resonance. Our conclusions are supported with semi- 
classical calculations of the relevant AGs splitting ruled 
by the PB resonance. 

This work was supported by the MIGINN- Spain un- 
der contracts iMath-COmOLIDER GDS2006-32 and 
MTM2009-14621, UBAGyT (X237) and GONIGET (Ar- 
gentina) . 



A final point worth addressing is the limit A^ ^ oo of 
our results, which is closely related with the extent to 
which quantum mechanics can distinguish the PB struc- 
ture of a given classical resonance. According to KAM 
and PB theorems such structures exist even for infinites- 
imal values of the perturbation, although as k decreases 
they span smaller regions of phase space. On the other 
hand, the finiteness of h imposes a limit to the phase 
space details which can be resolved quantum mechani- 
cally. How do they two extremes reconcile with our pre- 
vious conclusions is not clear in principle. For example, 
in the correlation diagram of Fig. [T]the 1:6 resonance is 
only observed for k > 0.166, point at which the first ele- 
ment, i.e. niower = 0, of the corresponding family of AGs 
is located. Gould this resonance be resolved for smaller 
perturbations if a smaller value of h was chosen? Fortu- 
nately, this question can be easily answered in our cal- 
culations, where is a parameter directly related to the 
size of the Hilbert space, N = {27Th)-^. Indeed, the an- 
swer to the question is yes, as demonstrated by Fig. [3l 
where the results for for three different values of 
are shown. For increasingly larger values of N the first 
member of the AGs series appears at smaller values of /c, 
the interaction being satisfactorily described in all cases 
by the semiclassical theory. 

Finally, let us remark that our conclusions are also 
valid for other resonances and systems. This is illustrated 
in the inset to Fig. [3l where the results for the resonance 
1:10 are shown, and a perfect agreement between quan- 
tum and semiclassical results is also seen. Moreover, a 
similar study has been carried out in a realistic model for 
the vibrational dynamics of the LiGN molecule [23] arriv- 
ing at the same conclusions; the results will be reported 
elsewhere. 

Summarizing, in this Letter we have shown that quan- 
tum manifestations of the PB theorem, a cornerstone to 
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